MAXIMUM MODULUS ESTIMATE FOR THE SOLUTION OF THE 
NONSTATIONARY STOKES EQUATIONS 

TONGKEUN CHANG AND HI JUN CHOE 

CN . Abstract. A maximum modulus estimate for the nonstationary Stokes equations 

in C^ domain is found. The singular part and regular part of Poisson kernel are 
analyzed. The singular part consists of the gradient of single layer potential and the 
gradient of composite potential defined on only normal component of the bound- 
j^ , ary data. Furthermore, the normal velocity near the boundary is bounded if the 

boundary data is bounded. If the normal component of the boundary data is Dini- 
continuous and the tangential component of the boundary data is bounded, then 

0^ ' the maximum modulus of velocity is bounded in whole domain. 

(N 



O 



(N 

in 



Keywords: Nonstationary Stokes equations, Maximum modulus, bounded cylinders. 



< 
c^ ■ 1. Introduction 



A maximum modulus estimate of the nonstationary Stokes equations is presented. 
In the case of the stationary flow, Maremonti and Russo[l] obtained a quasi maximum 
principle and Varnhorn[TO] showed a maximum modulus theorem for C^'° domain: 

max|ti| < C(r2) max |m|, 

Q I where m is a solution to the stationary Stokes equations in domain Q. We also note 

^ . that Maz'ja and RossmannjS] considered the maximum modulus estimate for the 

stationary Navier-Stokes equations in polygonal domain. In a canonical domain like 
K> ' ball, Kratz[3] found the best constant C{Q) such that 

- - - max \u\ < —n[n + 1) max \u\, 

xeBi 2 xedBi 

where Bi is the unit ball in R". For a more general domain like Lipschitz in R^, 
Shen[7j obtained a maximum modulus estimate and the higher dimension problem is 
still unresolved. 

The maximum modulus estimate of the nonstationary problem is heavily entangled 
with the structural form of Poisson kernel and the solvability of the boundary value 
problem is essential. As a classical result, Solonnikov[8] solved the initial-boundary 
problem in C^ domain for the isotropic Sobolev spaces and later he [9] extended the 
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solvability to the anisotropic Sobolev spaces. The L^ solvability for the Lipschitz do- 
main was obtained by Shen[6] for any dimension and Choe and Kozono[T] considered 
the case for the mixed norm potential spaces. 

To be more specific, we state the nonstationary Stokes equations: 
Ut-uAu + Vp = infix(0,r), 

divu = infix(0,r), 

(1-1) , 

u\t=o = in f2, 

u\dnxio,T) = g ondnx (0, T), 

where Q is C^ bounded connected domain in R" and < T < oo and u is the 

viscosity which we assume 1. In addition, we assume the boundary data g satisfies 

the compatibility condition: 

g . Nda = 



'an 
for almost all t, where A^ is the outward unit normal vector on the boundary. Since 

nontrivial initial data can be treated by solving homogeneous boundary value prob- 
lems, we consider only the initial-boundary value problems with zero initial data. 

Contrary to the stationary case, the quasi maximum principle fails, namely, there 
is an unbounded solution whose boundary data is bounded. Heuristically speaking, 
at the boundary point where the normal component of boundary data has a jump 
discontinuity along an [n — 2)-dimensional surface on the boundary passing to it, the 
tangential component of the velocity blows up in the neighborhood of it. So we can 
not expect the quasi maximum modulus theorem like the stationary case. 

In this paper, we only consider the case that the space dimension is greater than or 

equal to 3. Dimension 2 case follows exactly the same path with logarithmic kernels. 

Denote E for the fundamental solution to Laplace equation and F for the fundamental 

solution to heat equation with unit conductivity. For a given boundary point y G dQ 

N[y) is the outward unit normal vector at y. We define the (n — 1) -dimensional 

convolution 

S(/)(x)= [ E{x-y)f{y)da{y) 
Jan 

for real- value function / : R" — )■ R which is just the single layer potential of / on 
dQ. We need a composite kernel. We define a composite kernel function k{x, t) on 
n X (0, T) by 



f dr 

i^{^^t)= r.^j, A x-y,t)E{y)da{y) 
Jan dN{y) 
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and a surface potential T for / by 

T(/)(x,t)=4/ / K{x-y,t-s)f{y,s)da{y)ds, 
Jo Jan 

for real-value function / : R""*^^ — )■ R. We state our main theorem: For given x G f2, 
X is the nearest point of x on dVt such that dist{x,dQ) = \x — x\ and for a vector 
valued function v{x), we define the normal component and tangential component to 
the nearest point x by 

vn{x) = {v{x) ■ N{x))N{x) and vt{x) = v{x) — {v{x) ■ N{x))N{x). 

Theorem 1.1. Suppose that the domain Q is bounded C^ and u is a solution to U.l\) 
for bounded boundary data g. The normal component of the velocity un is bounded 
and there is also a constant C{Q) such that 

max \uNix,t)\ < C(Q) max \q(y,t)\. 
{x,t)enx{o,T) {y,t)£dnx{o,T) 

Furthermore, the tangential component of the velocity u satisfies that 

r max |MT(x,t)-VS(^-iV)T(x,t)-VT(^-iV)r(x,t)| <C(fi) max \giy,t)\ 

Define the modulus of continuity of / at xhj uj{f){r,x) = supj^g^^^^^p,,^ \f{y) — f{x)\ 
and we say / is Dini-continuous in Q if 



\Dini 



r° dr 

,n = sup / oj{f){r,x) — < oo 
x&n Jo '^ 



for an tq > 0. From a direct computation, we have VS(/) and VT(/) are bounded 
if / is Dini-continuous on dQ and we obtain a maximum modulus estimate: 

Corollary 1.2. Suppose that the domain Q is bounded C^ and u is a solution to 
( fi.ij) . Suppose g is bounded on dQ and the normal component g^ is Dini-continuous. 
Then, there is a constant C{Q) depending only on Q such that 

max \u{x,t)\ <C{Q){ max \g{y,t)\ + \\gN\\Dini,dn)- 
(x,t)enx(o,T) (j/,t)69nx{o,T) 

As a separate interest, we obtain an improved L^ theory like Lemma 4.1. When the 
L'^(dQ) norm of the boundary data is bounded in time, then ||m(-, t)||j^2(Q) is bounded 
in time. Consequently, the local boundedness holds too(see Corollary 4.2.). 
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2. Kernels on half plane 

To study the equation (11. ip . we consider the case oiQ = R" = {(x', x„) G R" | x' G 
R'"~^, < x„ < oo} and for the notational simphcity we set D^^ = -^ and double 
indices means summation up to n. For notation, we denote x = (x',x„), that is, 
x' = (xi,X2, ■ ■ -jXn-i)- Indeed, the symbol / means the coordinate up to n — 1 and 
Un is the surface area of the unit sphere in R*^. 

We let r be the fundamental solution to the heat equation such that 

{|:r|2 
^=Tre"^F t>0 
0, t <0 

and H be the Newtonian potential of F such that 

H{x,t)= f T{y,t)E{x-y)dy. 
The Stokes fundamental matrix (F, 7) for R", n > 3 is 

Fij{x, t) = 6i,r{x, t) + _^ D,^D,^H{x, t) 
6(t) Xi 



where 5{t) is the Dirac delta function and 5ij is the Kronecker delta function. 
The Green's matrix (G, 9) for the half space R" is 

Gij{x,y,t) = 5ij{T{x-y,t) -T{x-y*,t)) 

+ 4(1 - 6jn)D,^ f " / D,Mx - ^)r(^ - y*,t)dz 

Jo JR"-1 

6j{x,y,t) = {l-6jn) { / Dx^E{x' - z',Xn)T{z' -y',yn,t)dz' 

V^R"-i 

+ / E{x' - z', Xn)Dy,r{z' - y', yn, t)dz' ) , 

where we denote x* = (x', — x„). 

The Poisson kernel {K, vr) for the half space is defined by 

Kij{x' - y',Xn,t) = ^^|£^|j;„=o - SjnOi{x,y,t)\y„=o 

= -26ijD^^T{x' -y',Xn,t) + ALij{x' -y',Xn,t) 

-6jnSit)D.^E{x'-y',Xn), 
7i,{x' - y', xn, t) = -25{t)D.,^D^^E{x' - y', x^) + AD,^D.,^A{x' - y\ x„, t) 

+ADtD,A{x'-y',Xn,t), 



MAXIMUM MODULUS FOR STOKES EQUATIONS 

where we defined that 

L,j{x, t) = D,^ f f D,J{z, t)D,^E{x - z)dz, 

Jo JR"-1 

A{x, t)= T{z', 0, t)E{x' - z', Xn)dz'. 

Lij and A satisfy the estimates 
(2.1) \DtDlwrU,{x,t)\< 



(2-2) \DiDTA{x,t)\ < J , 

where 1 < i < n and I < j < n — 1 (see [2] and [8]). The estimates (12. ip of Lij and 
the estimate of the fundamental solution to heat equation F imply that 

(2.3) \DlD'jDrK,j{x,t)\< 



1„ I 1 ; 



The solution {u,p) of the Stokes system (11.11) in fi = IV^ with boundary data g is 
expressed by 

u'{x, t) = Y,]=i /o Inn-i Kij{x' - y', Xn, t - s)g-j{y', s)dy'ds, 

\ ) f 

p(x, t) = Yl^^^i Jo Inn-i TTjix' - y', Xn, t - s)gj{y', s)dy'ds. 
We have relations among L and A such that 

(2.5) / ^ La = —ZDx^i- , Lin = Lni + Bin, 

l<i<n 

where Bin{x,t) = J-^„_^ Drc^r{x' - y',Xn,t)Dy^E{y',0)dy' = ^^K{x,t) \i % ^ n and 

-Bnn = 0. 

For further computation, we introduce Lebesgue spaces and Sobolev spaces: 

L^(fi) = {/; / \j\Ux < oo}, W'^iQ) = {/; [ \f\^ + \Vf\^dx < oo}, 
Jn Jn 

L^(0, T- W^'^{n)) = {f; [ [ l/r + \Vf\Pdxdt < oo}. 

Jo Jn 

3. Maximum modulus estimate in the half space 

In this section, we consider the maximum modulus estimate in the half space. 
The normal derivative D^^T has uniformly bounded L^ norm with respect to Xn on 
9R" X (0,r)(see (13. 9p ) and hence we focus only on the kernel function Lij. By 
introducing a composite kernel n we are able to identify the singular kernels. The 
following lemma is a key stone for the maximum modulus estimate. 
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Lemma 3.1. Let 1 < i < n and 1 < j < n — 1. Then 

(3.1) / / \Lij{x',Xn,t)\dx'dt<C, 

Jo >/r"-i 

where C > is independent of Xn > and hence it follows that 

/•oo /• 

(3.2) / / \Lin{x',Xn,t) - Bin{x',Xn,t)\dx'dt < C, 
Jo JR"-! 

where C > is independent o/x„ > 0. 

The maximum modulus theorem for the half space follows from Lemma 13.11 

Theorem 3.2. Let g = {gi,g2,- ■ ■ , Qn) G L°°((9R" x (0, T)) and {u,p) is represented 
by (I23D. Then, 

(3.3) \\ut — VST(fi'n) — VTr(5r„)||^oo(Rnx{0,T)) < C'||fi'||L°°(9R!J:x{0,T)) 

for some C > 0. Furthermore, the normal component of the velocity u is hounded and 
there is also a constant C such that 

max \un{x,t)\<C max \q{y,t)\. 

(a:,t)6R!fx(0,T) (j/,t)eOR!;:x{0,T) 

To show the L^ boundedness of Lj^, we note that 
(3.4) 

Lij{x, t) = 2^n^ / t'^y^e"^ / DyT'{y', t)Dy^E{x' - y', x„ - yn)dy'dyn, 

Jo 7R"-i 

where F' is Gaussian kernel in R"^"^. 

Lemma 3.3. For 1 < j < n — 1, we get 

(3.5) 
I Iw-y'\<^W\ Dyr{y',t)Dy,E{^' - y\yn)dy'\ < Cr^e-^lx'!"' + Crt-i|x'|e-^ 

li'|l-'l<|j/'l<2k'|,|x'-y'|>i|x'|^%r'(y',t)Dj,„^(x'-2/',2/„)d2/'| <crt-5|x'|e- 

I k'\<\\A Dyr{y\t)DyMx' - y\yn)dy'\ < C|xr"4,|<.^ \yfe^^''^"dy' 

I 4,|>,|^,| Dyr{y',t)DyM^' - y',yn)dy'\ < Ct'^ /^^|^,| |y'|-"+'e-l^'l^d2/', 

where C > is independent of x', yn and t. 



.kl^ 
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Proof. Using integration by parts, we get 



/ DyT\y', t)D^^E{x - y', yn)dy' 

JW-y'\<h^'\ 

(3.6) = / ^lZyi^T'iy',t)Dy„Ei^' -y',yr.)aidy') 

J\x'-y'\ = ^\x'\ \X — y \ 

- / ^'{y\t)Dy^Dy^E{x -y\yn)dy. 

J\x'-v'\<h\x'\ 



For y' with |x' — y'l = ||x'|, we get |r'(y', t)| < Ct 2 e at aji(\\Dy^E{x' — y' ,yn)\ < 
C ^ — TT^rr. Here, the first term of the right hand side in (13.61) is dominated by 

\x'\- 



L,_,. iu,| \T\y'MDyMx'-y\yr:)\a{dy') < Cr — e~— ——^ 
(3.7) ' 



(I^T+?/^)" 



I /i2 
n-1 \x'\^ , /ii 



<Ct- — e-^\x\ 

Since T ,_ n^i,,.Dy.Dx^E{x' — y',yn)dy' = 0, using the Mean value theorem, the 
second term of the right hand side of (13. 6p satisfies 

iix-y|<i|x'|(r'(y'>^) - T'{x',t))Dy^D,^E{x' - y',yn)dy' 
(3.8) - C\x'\t-—-^e-— ij,,_y|<i|,,| (|,,4/'7+i)'^+^'^^' 

By dM]) - ([31]), we obtain ( 133]) ^. 

For ( l3.5P o. note that for y' satisfying \\x'\ < \y'\ < 2\x'\ we have \x' — y'\ > ^\x'\. 
We have \DyT'{y',t)\ < Cr^^-il^'le"^ and Dy^E{x' - y',yn) < C|x'|-V, and 
thus we get 

f DyT'{y',t)Dy„E{x' -y',y^)dy' <Cr^~-2\x'\e-^. 

^\x'\<\y'\<2\x'\,\x'-y'\>^\x'\ 

Hence, we obtain ( l3.5P o. 

Since J. /|<i|^/| Dy.T'{y',t)dy = 0, using Mean-value Theorem , (l3.5P n is proved by 

/ DyT'iy', t) (d.,^E{x' - y', ?/„) - D^M^' , yS)dy' 

J\y'\<^\x'\ ^ ^ 

/|2 , \-il f +-^l /|2 --1^ 



J\y'\<\W\ 

<c(|xf + 2/„)-^ / , \y'?e-\y\'dy'. 



\y'\^-.^ 
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Finally, (13. 5p ^ follows by 

/ DyT\y',t)Dy„Eix' - y',yr.)dy' < Cr^ [ \y'r^'e-'-^dy' 

J\y'\>2\x'\ J2\x'\<\y'\ 

= Ct-"^ [ \yT^'e-\''^'dy'. 

n 

Following a similar proof to Lemma 3.3, we get the following lemma. 
Lemma 3.4. For 1 < i, j < n — 1, we get 

[ DyT'{y',t)Dy^E{x' - y',yn)dy' < Crt-i|x'|e-^ 

f DyT'{y',t)DyMx' - y',yn)dy' < Cr'^~-^\x'\e-^ , 

J l\x'\<\y'\<2\x'\,\x'~y'\>^\x'\ 

f DyV'{y\t)DyMx - v\Vn)dy' < C\xT [ , „ Ivfe-^y'^'dy' , 



f DyT'iy', t)DyM^' - y', yn)dy' < cr"^ [ , 

J\y'\>2\x'\ J^<\y'\ 



\y'\-'''^^e-\y'\ dy'. 

Proof of Lemma 13.11 

Note that for 1 <i <n and 1 < J < "«- — 1 

/ / \Kij{x ,Xn-,t)\dxdt < I I \Dx„T{x' ^Xn-,t)\dx'dt 
Jo Jr"-! io Jr"^ 

+ / \Lij{x',Xn,t)\dx'dt. 

Jo Jr"-i 

2 

Here, using change of variables (^ = s), we get 

/o /r"-i l^^nr(x', Xn, t)\dx'dt = C/g t-2x„e"^ j^„_^ t-'^e'^dx'dt 
(3-9) = Cxr. J^ t-le-'f J^„_, e-\-'\'dx'dt 

Hence, to prove Lemma [3?T| it is sufficient to show J^ j-Rn-i \Lij{x\ Xn, t)\dx'dt < cxo 
for 1 < z < n, 1 < j < n — 1. 

By (12. ip ^. for 1 < z < n and 1 < j < n — 1, we get 

/;' /j^„_, \Lij{x',Xn,t)\dx'dt < Cjf- J^„_, t-5(|xf + Xl + ty^dx'dt 

<cjf-H{xl + t)-Ut = c. 
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To calculate J^j /rk-i \Lij{x' ,Xn,t)\dx'dt, we may assume x"^ < T. By the represen- 
tation (13 ■4p . and Lemma [3.31 and Lemma [3 ■4[ we have 

(3.11) ^ 

Jx"^ Jj^ri— 1 \''-'ij\-^ , Xn,l^ )\CIX UZ 



-^Jx2jRn-iJo ^ '?/ne ' [^ ' e ' |X I ' + t 2 2|x |e t 



-2 yt Vt 

I + 11 + III + IV, 



where 



T r rx 



3 _H^ ra-1 _|a:'|2 



1 = / ^ ^2^"^ * ^ ^ ^ ' P I di/ndx at, 

Jxl JR"-1 JO 

11= / / ^~^y„e ?'t~^~2 |a;'|e ~dyndx'dt, 

Jxl JR"-! io 

3 a4 



^2 jRn-l Jo 'f\y'\<^^ 

IV= r I r t-ly^e-'ft-'^ [ \y'\-^'e-\y'\'dy'dyMdt. 

Jxl jR"-i Jo 7aM<|y'| 

Using change of variables twice, we have 

I = /^2 t-^~^ /r„-i e-^|x'|-^t/o^ ynt-y^dyndx'dt 



(3.12) ^ C/J^t-5-VK„-.e-^|x'|-H(^)Wdt 

= C 



and 



// = j'^2t 2 Vj^„_ie "t |x|t/o^?/„e y^^dyndx'dt 



,'|2 



(3.13) < C/J^ri-VR.-.e-^|x'|t(^)Wdt 
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We divide /// into two parts /// = I Hi + III2, where 



1 — 


r 


1 


rXn 2 

/ ,_3 _an ,, 
/ t 2y^e t X 


— n 


/ ,., 




Jxl . 


'|x'|<v^. 


Jo 


- 


Vi<l^ 


2 = 


f 


/ 


/ t 2y„e t x\ 


— n 


/ ,. 




Jxl. 


'|a;'|>v^. 


/o 


- 


Vl<4^ 



ly'l^e '^ ' dy'di/ndx'dt, 
\y'\^e~^y'^'dy'dyndx'dt. 



Here, 



II h < c 


rt-l [ x' -(^)"+i r y^e-'fdyMdt 

Ixl J\x'\<Vi Vt Jo 






l-T f 


= c 


/ xlr-i-^ / x'\dx'dt 




'xl J\x'\<Vt 




fT 


<c 


/ xlr^dt 

/t2 


= c, 


^n 


II I2 < c 


/ t"^ / x' ~" / yne~~t^dyndxdt 




'xl J\x'\>y/t Jo 




r^ . r 


<c 


/ t-2x2 / x'\-'^dx'dt 




Ixl J\x'\>Vt 




fT 


<c 


I r^xldt 



= c. 

Hence, we get 

(3.14) /// < C. 

Similarly, we divide IV into two parts 



JxlJ\x'\<^J0 J xl J Vt<\x'\ J 

= IVi + IV2. 
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Here, with straightforward integrations 

rT /• f-Xn 2 

IV, <C t-t-2 / / y^e-'-^dyMdt 

Jxl J\x'\<VtJo 

<C t'^'^xl / dx'dt 

Jxl J\x'\<Vt 

<c [ r^xldt 
= c 



and 



IV, <C t-t-f / / y^e-'f / \yr''^\-\y'\'dy'dyMdt 

Jxl Jvt<\x'\ Jo ■^^<\y'\ 

< C Tt-t-fx^ f [ ly'l-'^'-'e-^y'^'dy'dx'dt 

Jxl JVt<\x'\ J'-^<\y'\ 

= C f r'^xl f f \y'\-^+\~\y'\"dy' dx'dt 

Jxl Jl<\x'\ J\x'\<\v'\ 



<c. 

Hence, we get 

(3.15) IV <C. 
Therefore, from (13.101) - (13.151) . we prove 

(3.16) / / \Lij{x',Xn,t)\dx'dt<C 



JR"-1 

for 1 < 2 < n and 1 < j < n — 1, where C is independent of Xn- With (13. 9p . this 
imphes (13. ip . 

By the second identity of (12. 5p and (I3.16p . we prove (13. 2p for the case i ^ n, and 
by the first identity of (12. 5p and (I3.16p . we prove (13. 2p for the case i = n. This ends 
the proof of Lemma 3.1. 
D 
Proof of Theorem 3.2. 
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We begin the proof of Theorem 3.2 by the representation (12. 4p of u such that 
u\x, t) = y2 / ^ij(^' ~ y'^ ^n, t - s)gj{y', s)dy'ds, 

n-l 



n-i „t p 

y2 Kij{x' -y',Xn,t- s)gj{y',s)dy'ds 

~[ Jo JR"-1 

/ / Kin{x' -y',Xn,t- s)gn{y',s)dy'ds 

Jo JR"-1 



and the last potential for Qn is written as 



/ / Kin{x' -y',Xn,t- s)gn{y',s)dy'ds 

Jo JR"-! 

= - 26in / / D^^r{x' - y\ Xn, t - s)gn{y', s)dy'ds 

Jo JR"-i 

+ 4/ / Lin{x' -y',Xn,t- s)gn{y',s)dy'ds 

Jo JR"-i 

- -— / / E{x' -y',Xn)gn{y',s)dy'ds. 

(y^i Jo Jr."-^ 



Since Lj„ = L„j + Bin, we have 



/ / Lin{x' -y',Xn,t- s)gn{y',s)dy'ds 
Jo 7r"-i 

Lniix' - y', Xn, t - s)gniy', s)dy'ds 
■^— K,{x' -y',Xn,t-s)gn{y',s)dy'ds 

OXi Jo jR,n-l 



JR"-i 

d r 



for 1 < i < n— 1 , where we defined the composite kernel function k{x, t) on R" x (0, T) 
by 

f 9 

K{x,t) = / — — T{x' — z' ,Xn,t)E{z' ,0)dz'. 

Define the surface potential T{gn) by 
(3.17) T{gn){x,t) = A [ [ k{x' - y' ,Xn,t - s)gn{y' ,s)dy'ds. 

Jo JR"-i 

Moreover, we have that 

-— / E{x' - y\ Xn)gn{.y\ s)dy' = ^— S(fi^„) 

(y^i JR"-1 OXi 

Therefore we conclude that the tangential part, which is associated with Lij, satisfies 

\uT{x,t) -VSTign){x,t) -VTT{gn){x,t)\ < C||fi^||Lcx>(R,n-ix(o,T)) 
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for all (x,t) G R!^ X (0,T). 

The normal velocity w„ behaves even better. First, we know that -^S{gn) is the 
Poisson kernel expression of the solution for the Laplace equation in the half space 
and satisfies the maximum principle. In the case i = n, we have a relation from (12. 5p 

l<i<n-l 

which has a bounded V" norm on the lateral surface. This conclude the maximum 
modulus estimate of Un- 

U 

4. Maximum Modulus Estimate in C^ Domain 

We denote the Green's matrix for the domain Vt by (G*^, 9^) and for a given point 
X G f2 we let x G dVt satisfy |a; — a;| = dist{x,dQ). The interior L°° bound estimate 
can be shown by the layer potential method in [6] and we consider separately the case 
that the generic point x is close enough to dQ. 

Indeed, to see the interior boundedness, we need to show the boundedness of the 
double layer potential in L°°{0, T; L'^{dQ)). Since the boundary data is bounded, we 
can represent the solution by the double layer potential in [6] from L^ theory such 
that 

(4.1) Ui{x,t) = / '/ {x-y,t- s)hj{y, s)da{y)ds 



an dNiy) 



--{Kh)i{x,t) 



:h{y,t)-N{y)da{y) 



and 



(4.2) g = J_h + ^h={-h + K)h 

in the sense of L'^{dn x (0,T)) for an /i G L'^{dVt x (0,r))(see Theorem 2.3.6 and 
Theorem 5.1.2 in [B]). Furthermore — |l + K is invertible on L'^^{dVL x (0, T)), where 
the subscript a means solenoidal. From the representation, we have a continuity 
lemma in time of the density function h in (14. ip . 

Lemma 4.1. The inverse of the double layer potential — |l + K is hounded in time 
as an operator of L'^{dQ) and there is a constant 6 > such that if \ti — t2\ < 6, 

\\hi-,h)\\LHdn) < C||fi'(-,ti) - gi-,t2)\\L^dn) + C\\h\\L^^^ofy;L2{n))- 
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By an iteration there is C such that 

ll(~2^~'' "^'^ 5'l|L°°{0,T;L2(n)) < C\\g\\L^{Q^T;L-2(n))- 

Proof. We assume the boundary data g G L°°(0,T; L?'{dVL)) and after arranging the 
singular integrals in the double layer potential expression we have 

gi{x, t2)-gi{x, ti) = --{hi{x, ^2) - hi{x{t - 1)) 

^■'""'^ -ih{y,t^)-h{y,t2))■N{y)da{y) 



an ujn\y-x\ 



*2 r Qp. . 

I dm ) ^^~y'^^^ s)hj{y, s)da{y)ds 

= {-h + H)(/i(-, t2) - h{-M)) + Ei/i + Es/i 

for almost all < ^2 < ^2 < T and x G dVt. 

We claim — |l + H : L'^{dQ) — )• L'^{dVt) is invertible and 

||(--I + H)"^e||i2(9f^) < C||e||i2(9n) 

for a constant C. First of all, if we set e = (— |l + H)/ and consider the normal 
components, then we have 

where A^ ■ H is the standard double layer potential operator of Laplace equation and 
— |l + A^ ■ H is invertible. So given vector valued function e G L'^{dQ), there is a 
scalar valued function w G L'^{Q) satisfying 

with H'wlliaj-gf^) < C||e||L2(gf^). Here w is the normal component of / and the tangen- 
tial component v of / is obtained by 

V = -2(e - cnN) - 2(iiw - {N ■ U)wN). 

Therefore we get 

f = v + wN 

and / satisfies 
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It remains the estimate Ei/i and £2^. Since Q is C^ domain, in the case of Gaussian 
kernel, there is C such that for all (x, y, t) G dQ x dQ x (0, T) 

d \x — up i^-yp 



Therefore we get from Minkowski inequality and Young's convolution inequality 



an 



dV 



an 



dNiy) 



[x -y,t2- s)h{y, s)da{y)ds 



da{x) 



<C 



<C 



< c 




ti \/^2 — S yjt^ — 5" ]qq_ ti — S 



\x-y\ 



__ I -^^ ^ I 



do{x) 




\x — y\ 



-e 2(t2-»)|/i(y, s)|rfo-(y) 



daix) 



\\K--.s)\\L-^{Si)ds 



< C\/t2 — tx\\h\\i^^(t^^tr,L'i{a))- 
By the same token, assuming ||/i||Loo(o,ti;L2(n)) is bounded, we have that 

||E2/l(-,tl)||L2(gQ) < C||/?.||L°°(0,ii;L2(f7)). 



D 



We let the generic point x be away from the boundary, say dist{x, dQ) = tq > 0. 

Since the kernel of the double layer is bounded by „5i_, for each e > and the 

''0 

density function h oi g for the double layer potential is bounded in L°°(0, T; L'^{dfl)), 
the interior L°° estimate follows. 

Corollary 4.2. Suppose the boundary data g G L°°(0,T; L^(9r2)). If dist{x,dVt) > 
ro > 0, e > and t < T, then there is C such that 

I / \ I ^ 1 1 1 1 

\u[X,t)\ < „_i_J |g'||L°°(0,T;L2(8n))- 

Now we start the boundary estimate. Since Stokes equations is translation and 
rotation invariant, we assume that x = and x = (0,a:„), x„ > 0. If x is close 
enough to dQ, there is a ball Bj.{0) centered at origin and C^ function $ : R"^^ — t- R 
such that ^ n Br{0) = {x„ > ^{x')} fl Br{0). Furthermore, $ satisfies that 



(4.3) 



$(x')| < C\x'\\ |V'<l>(x')| < C\x'\, |VV$(x')| < C 
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for x' G -Br(O) and the outward unit normal vector A^(a;', $(x')) at (a;',$(a:')) G 

af]n5,.(o) is 

iV(x',<l>(x')) = ^ =(y'^[x'),-\). 

We define a transform /i : f] fl i?r.(0) — t- R" such that 

Ky) = f^iy', yn) = {y, yn - ^{y')) 

and note that ii^^{y',yn) = {y',yn + $(?/'))• Since our generic point x is (0,x„), we 
have iJ,{x) = x. Hence the Green's matrix G on the half space can be transformed to 
a function fiG on Q such that 

fiG{x, y, t) = G{fi{x),fi{y),t) = G{x', x„ - <l>{x'),y', y„, - <!>{y'),t) 

and satisfies the zero boundary condition 

^IG{x,y',^y'),t) = 0. 

Moreover, the transformed Green's matrix {fiG, fi6) satisfies a perturbed Stokes equa- 
tions in fi X (0,T) 

d d 

— (/iG')y - Ay{nG)ij + -Q-if^0)i 

- Dy„Dy,X^^G).JDy^^Dy^<^ - D y^ (flO) , D y ^ ^ 

= 6ijd{x - y)d{t) + R{x, y, t) 

and the solenoidal condition 

Dy^{fiG)ij = -Dy^{fiG)ijDy^^ = Si{x,y,t). 

Therefore, if we let the perturbation (J, 77) = {G^—fiG, 6^—fi6), then (J, 77) satisfies 
the perturbation equations: 

d d 

(4.4) —Jij{x,y,t)-AyJij{x,y,t) + —rjj{x,y,t) = Rij{x,y,t) 

(4.5) Dy^Jij{x,y,t) = Si{x,y,t), 
where R is 

i?,, = - D,„(/xG'),,A;$ - 2D,,D,„(^G),,D,,$ 

+ Dy„Dy^itiG)ijDy^^Dy^^ + D,„ ( ^^ ) , Z^,^ ^ 

=1 + 11 + III + IV. 
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We have already discussed the boundedness of velocity u in the interior by double 
layer potential in L^ theory, we begin to prove the boundedness near the boundary. 

The plan to get L^ bound of the perturbation J of Poisson kernel on dQ x (0, T) 
relies on the Lp(0,T; VF^'^(r2)) estimate and the trace theorem for it. Recall that 
the Poisson kernel is a derivative of Green's matrix and that is the reason that we 
need ^^(0, T; W'^'^{Q)) Sobolev type estimate. Therefore, we need to estimate the L^ 
norm of R in (jO]), Lp{0, T; W^'P{n)) norm of S and Lp{0, T; W-^'P{n)) norm of St 
in (fii nft) X (0, T) in (i3]), where W-^^p{Q) is the dual space of W^^p{Q). 

Since the Green's matrix G is associated with the Gaussian kernel and the com- 
posite kernel H, we estimate their derivatives first. We have 

\DyJ{x -y,t)\< -^^JL_Ae ^ 

\Dy„Tix-y*,t)\ < —^ ^ g-JL^42i±£ 



2 

It 






\Dy,DyJ{x-y\t)\< r^ ^2 ^ 
Since \Dy^^{y')\ < C\y'l |A'<l>(y')| < C and x' = 0, we get 

\DyJix-y,t)A'^y')\ < C\y -x\ ^_ ^y'i^.^y.-..^^ ^ L^((f] n B.) x (0,T)) 

\Dy,DyJ{x-y,t)V'Hy')\ < ^+i '^ ~ f e" '"'''"""'^"'' e 1^(1] n 5.) x (0,r)) 

as a function of y for all p G [1, ^^)- In the same way, we have 

DyJix-y*,t)A'^y'), Dy^DyJix - y,t)V^y') E LP{{nnBr) x (0,r)) 

as a function of y for all p E [1, f^)- 
Applying (12. ip . we have 

\DyA. / / D^^E{x - z)V{z - y\t)dz\ < -^— ,^ 



Dy,Dy,D,^ I " / D,,i?(x - ^)r(z - y\ t)dz\ < - 

Jo JR"-i t2 



t^(|x'-yf + |2/„ + x„|2 + t)5(-+i)' 
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Hence, we have for p G [1, ^^) 

\Dy,.D,, r [ D^^x - z)T{z - y*,t)dzA'^y')\ 

Jo JR"-i 



< ^^1 , \, ^ ^ L^mnBr) X (0,T)) 



\Dy,Dy„D,^ [ " / D^M^ - z)T{z - y*,t)dzV^y')\ 

Jo JR"-i 



< -^^=2==^ e L^m n Br) X (0, T)). 

Although there is a transformation yu of domain, these estimates imply that /, // 
and III are in LP((fi n Br) x (0, T)) as a function of ?/ for all p G [1, ^). 

It remains to get L^ estimate of the pressure kernel 6. For each fixed time t, we 
have 

\IV\ <C\y'\\ [ D,,i?(x'-z',x0^4^e-^^4l^c/^'| 



(4.6) + C\y'\ 



D,, D,, / — n — i=^e 2t dz 

^R"-i -y/|a;' - 2:'|2 + x^ vc 



The first term on the right is U for p G [1, ^^) by the Young's convolution inequality 
since the kernel , ^" tt has bounded L^(R"~^) estimate as a function of z' inde- 
pendent of a;„. For the second term, we recall the following proposition by Solonnikov 
(Proposition 2.3 in [9]): 

Lemma 4.3. Let M(x, t) he a function defined for x G R" and t > and having the 
properties 

M{Xx, \H) = A™M(x, t), A > 0, 

\D',DtMix,t)\ < Ct^'^^exp U-^ 

Then the integral 

J{x,yn,t) = E{y)M{x' -y',Xn,t)dy' 

JR" 1 

satisfies the conditions 



JiXx,Xyn,XH) = X"'+^J{x,y^,t), 

m + n-l-2s-k \k'\+l + n-2 ^2 

\D'^DlDtJ{x,yr.,t)\ < Ct ^ " (|^'|2 + (a;^ + ^j2^^^- . ^--^ 

where k = {ki, ..., kn) and \k'\ = A;i + ■ ■ ■ + /c„_i. 
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See Proposition 2.3 in [9]. 

So, we find the second term of (14. 6 p witliout tlie transformation ^ is bounded by 

VtVW ^ Vt ' ^ J ^ 



,c(^,,)^(,.,,,.^,,,)--.-. 



n-l 



wliicli is in LP{{Vt n B^) x (0, T)) for p G (1, ^). 

Hi Ic In TP for all n C: H 

-' n+1 



Tliis concludes tliat 11/ is in U' and i? in (14. 4p is in L^ for all p G (1, ^^^) after 



adjustment of the domain transformation /i. 

To get L^(0, T; iy^'^(r2 fl Sr)) bound of S in (14.51) we follow a similar program to 
R. Indeed, we have 

VyS = -VyDy„if,G).,Dy^^y') - Dy,XfiG),,VyDy^^y'). 

The terms in the right hand side have already been considered in the estimates of J, // 
and /// of i? except Dy^Dy^{fiG)ijDy^<^{y'). But, Dy^Dy^{fiG)ij = T.i<j<n~iDy,Dy^{fiGij] 
and hence Dy^Dy^{^G)ijDy^(^{y') has the form of //. Therefore, we get 

\\S\\Lv{G,T;W^^v{nnBr)) < C independent of x. 

It remains to find Lp{0,T; W~^'^{Q fl B^)) estimate of DfS. Since S is defined as 

S,{x,y,t) = -Dy^Sf^G),,Dy^<l>{y') 

and $ is independent of y„, 1/^(0, T; W~^'^{Q fl 5,.)) norm of St is bounded by 

G [ [ \DtifiG)V^y')\Pdydt 
Jo JnnBr 

for a constant C. By disregarding $, we have 

AG,, = S,,{DtT{x -y,t)- DtT{x - y\t)) 

+ 4(1 - 5,„,) A^ / " / D^M^ - z)DtTiz - y\ t)dz 

Jo JR"-1 

and 



2 \ , ,2 /I U „,*|2 



1 \x — y\ \ \x-v\'^ ( 1 \x — y 



* "■' - ^ vT-^^ vr^v I vr^^ v/r^' 



\DtGi,\ <G\ ^^x?+ -„:. I e ~ + C I ^;:xo + _1V I e" ^^ 
+ C 



A, / " / A,^(x - ^) Ar(z - y*, t)dz 

Jo 7R"-i 
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and Proposition 2.5 in [9], we have 



D,^ / / D.,^E{x - z)DtT{z - y*,t)dzV^iy') 

Jo JR"-! 



C\y'\ vl 

< — we~^. 



This imphes that DtGijV^{y') E LP{{n n Br) x (0,T)) for all p E [1, ^). 

Since the estimates of R and S hold only in a small ball near boundary, we need 
a localization. For the localization, we take a cut off function (p such that = 1 in 
Br and = in the complement of B2r and we consider {(j)J, (prj) as a solution to 
the inhomogeneous Stokes equations. We delete the generic point x in the various 
expressions. Therefore by Theorem 3.1 in [S], we obtain the following lemma for the 
perturbation {J,ri). 

Lemma 4.4. There is a constant C depending on r and Q such that 

\\J\\LP{0,T;W'^'P{nnBr)) + \\v\\LP{0,T;W^'P{nnBr)) 

< C(l + IIG \\LP{0,T;W^-P{nn(B2r\Br))) + 11^ ||LP{Qn(B2r\Br)x{0,T))) 

/ora//pG(l,gf). 

By the trace theorem in W^'P{Q fl Br), the following lemma also holds. 
Lemma 4.5. There is a constant C depending on r, Q such that 

" ^^ Lp{o,T]W^~p'^{dnnBr)) ^^ Lp{o,T;W^~p'P{dnnBr) 

<C{1 + \\G \\LP{0,T;W^'P{nn{B2r\Br))) ~^ W^ \\LP{nn{B2r\Br)x{0,T))) 

forallpE{l,^^). 

The generic point x is in Br and hence the Green's matrix {G^,9^) has no singu- 
larity in the complement of Br as a function of {y, t). Therefore we have that for all 
pE [l,oo] 

lie \\LP{0,T;W^-P{nn(B2r\Br))) + 11^ \\LP[nniB2r\Br)xiO,T)) < C 

for a constant C depending only on p, r and ^2. 

Now we prove our main theorem. The Poisson kernel K^{x, y, t) satisfies 

^""(^'2/'^) = Tr^G''{x,y,t)-e''{x,y,t)N{y), for all {x,y,t) E fix9fix(0,T). 
dN{y) 

We have that 

G^ = fiG + J, e^ = fie + r] 
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and from Lemma 4.3 we know that V J and r] have bounded L^{{dQ fl Br) x (0, T)) 
norms independent of x since L*'(0, T; W "p'^i^dQ fl Br)) for p G (1, ^^) is embedded 
in L^{{dilr\Br) x (0, T)). So we need to consider only gj^, . fiG{x, y, t) and fi9{x, y, t). 
The L^{{dQr\Br) x (0, T)) bound of fi9{x, y, t) as a function of y' for the generic point 
X = (0, x„) can be obtained by Lemma 4.1 after considering coordinate transform //. 
From the definition of the transformation of /i and the local representation of the 
boundary dQ, we have that for y = {y', ?/„) = {y', $(?/')) ^ ^^ ^ ^r 

1 d 

v/l+ |V'$(2/')I^2/" 

^ V.,G{x,y\0,t)-V',^y') 



1 9 
:t— G(x, 2/', 0, t). 

Furthermore, we have already proved that the Li{{dQ fl Br) x (0,T)) norm estimate 
of the first term 



T 

JdnnBr 



^ V'G{x,y',0,t)-V'^y'] 



dy'dt < C 



for some C independent of x since |V'$(|/')| < c\y'\. 

By the expression of Poisson kernel K we have 

d 
--— Gijlx, y\ 0, t) = Kij{x' - y', a;„, t) + SjnViix' - y', 0, t) 
oyn 

= -2SijD^„T{x' - y', x„, t) + 4(Lij(x' - y', x„,, t) - 5jnBin{x' - y' , x„, t)) 

+ AbjnBin{x - y', Xn, t) - 6jnS{t)Dx^E{x' - y' , Xn) 

We know already that —2SijD.j.^T{x' — y',Xn,t) + 4(Ljj(x' — y',Xn,t) — 6jnBin{x' — 
y',Xn,t)) has L^ bounded norm as a function of {y',t). Therefore in the solution 
expression for of u we can write 

/ / SjnS{t - s)D^^E{x' - y', Xn)gj{y, s)day 

Jo JdnnBr 

d f 

O^i JdnnBr 

4 / / 5jnBin{x' - y\ Xn, t - s)gj{y, s)day 

Jo JannB,. 



'dnnB, 
and 



22 TONGKEUN CHANG AND HI JUN CHOE 

= —T{g„){x,t). 

If we denote e„ = (0, 1) which is the coordinate vector for j/„, we have that the 
component of boundary data Qn is 

gn = g-N{y)+g-{en-N{y)) for y e dQ n Br 

where g = {gi,g2,- ■ -ign) is the boundary data and hence we have 



V,. / E{x' - y\ Xn)gniy, t)day = \/,S{g ■ NX' 



+V, / {E{x' - y', Xn) - E{x' - y', x^ - <l>{y')))g{y, t) ■ (e, - N{y))d(Ty 

+V:r / E{x' - y', Xn)g{y, t) ■ (e„ - N{y))day, 

JdnnBr 

where S is the single layer potential operator and X is the characteristic function. 
Since x' = and $(?/') < C|y'p, 

/ \V,{E{x' - y', x^) - E{x' - y', x^ - Hy'M day < C. 

Then, by observing that 

\en-N{y)\<C\y'\ 
we find that VxE{x' — y', x^) ■ (e„ — N{y)) has bounded L^ norm as a function of y' 
and we have that 



sup 

x&nnBr 



Vcc I E{x' - y', Xn)g{y, t) ■ (e„, - N{y))day 

annBr 



<C\\g\ 



Loo(anx(0,T))- 



Similarly we find that Vk(x' — y',Xn,t){en — N{y)) has a bounded L^ norm as a 
function of {y',t) and we have that 



sup 

{x,t)ennBrx{o,T) 



/ Vk{x' - y, Xn, t - s)g{y, s) ■ (e„ - N{y))dayds 

JannBr 



< C\\g\\L°°{anx{o,T))- 
With the interior L°° estimate, localization with the small balls B^ and the pre- 
ceding kernel estimates of L^ bound, we prove our main Theorem 1.1. 
For Corollary 1.2, we observe that 

|V.S(^-iVX^^^^OI<C^/ ^^b(2/',$(?/'),i)-iV(y',<^(?/')l^?/', 

JnnBr \y I 

where x = (0,x„). Similarly, we also have 

|V.T(<?-iVX^^^^OI<C r/ -^^\giy',^y'),t) ■ Niy',^y'))\dy'dt. 

Jo JnnBr \y 
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The boundedness follows from the Dini-continuity oi g ■ N. 

D 
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